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ABSTRACT 

It is shown that the McCall-Hahn theory of self-induced transparency in coherent optical 
pulse propagation can be identified with the complex sine-Gordon theory in the sharp line 
limit. We reformulate the theory in terms of the deformed gauged Wess-Zumino-Witten 
sigma model and address various new aspects of self-induced transparency. 
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Self- induced transparency (SIT), a phenomenon of anomalously low energy loss in coher- 
ent optical pulse propagation, was first discovered by McCall and Hahn|I[] and the integra- 
bility of the SIT equation was demonstrated by employing the inverse scattering method[0. 
When phase variation is ignored in the case for a symmetric frequency distribution g{Aw) of 
inhomogeneous broadening, McCall and Hahn have proved an area theorem for pulse prop- 
agation. In the sharp line limit where the frequency distribution is sharply peaked at the 
carrier frequency Wq such that g{Aw) = 6{w — Wq), the SIT equation reduces to the well- 
known sine-Gordon equation and the 27r area pulse becomes a 1-soliton of the sine-Gordon 
theory. However, when phase variation is included, the area theorem no longer holds and 
the structure of SIT in general has not been well understood except for the construction 
of explicit solutions by the inverse scattering method^]^. In particular, despite its inte- 
grability, the SIT theory in the sharp line limit has not been identified with a known 1-1-1 
dimensional integrable field theory, which made a systematic understanding of SIT in terms 
of a lagrangian field theory impossible. 

The purpose of this Letter is to show that the SIT theory with phase variation can be iden- 
tified with the complex sine-Gordon theory in the sharp line limit. The complex sine-Gordon 
theory, a generalization of the sine-Gordon theory with a phase degree of freedom, can be 
reformulated in terms of a nonlinear sigma model which is known as the integrably deformed 
gauged Wess-Zumino-Witten(WZW) model associated with the coset SU(2)/U(1)[|I|1. This 
allows us to address various new aspects of SIT in terms of characteristics of the complex 
sine-Gordon theory; e.g. topological v.s. non-topological solitons, local gauge symmetry, the 
U(l)-charge conservation, the chiral symmetry and the Krammers-Wannier duality for dark 
v.s. bright solitons. We also explain the off- resonance effect and inhomogeneous broadening 
of SIT in the context of the local gauge symmetry of the present formulation. 

The SIT equation is given by 

dE + 2p <P> = 
dD - E*P - EP* = 
dP + 2iAwP + 2ED = (1) 

where Aw = w — wq , d = d/dz , d = d/dz,z = t — x/c, z = x/c. E,P and D represent the 
electric field, the polarization and the population inversion respectively. The bracket denotes 
an averaging over the distribution function of inhomogeneous broadening. Since Eq.(l) is 
invariant under the interchange {f3,E,P,D) ^ {—i3,E,—P,—D), we assume the coupling 
constant (3 to be positive which we set to one by rescaling E, P and D. In a simpler case 
where phase variation is ignored to make E real and the frequency distribution is sharply 
peaked at the carrier frequency(Aw = 0), we may parametrize E, P and D by 

E = E* = dip , < P >= P = - sm2ip , D = cos2ip . (2) 

Then the SIT equation reduces to the well-known sine-Gordon equation 

ddip-2psm2ip = 0. (3) 
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In order to include phase variation as well as the off-resonance effect {Aw 7^ 0), we assume 
E to be complex and require that the distribution is sharply peaked not necessarily at the 
carrier frequency, i.e. g{Aw) = 6{Aw — ^) for some constant ^. Introduce a more general 
parameterization of E, P and D in terms of three scalar fields (f, 9 and 77, 

E = e'^^-^''\2dr]^^ - idip) , P = ie'^^-^'') sin2G9 , D = cos2ip . (4) 
sin 

The main result is that the SIT equation given in Eq.(l) then changes into a couple of second 
order nonlinear differential equations known as the complex sine-Gordon equation, 

— cos 09 — 

dd^ + A — :T^dr]dr] -2 sm2ip = (5) 
sin 

2 - 

ddn idridip + dridip) = (6) 

sm 2{p 

together with a couple of first order constraint equations, 

2 cos^ Lfdr] — sin^ LpdO — 2^ siv? = (7) 
2 cos^ i^dr] + sin^ i^dO = . (8) 

The complex sine-Gordon equation first appeared in 1976 in a description of relativistic 
vortices in a superfluid |0, and also independently in a treatment of 0(4) nonlinear sigma 
model0]. Note that Eqs.(5) and (6) consistently reduce to the sine-Gordon equation when 
phase variation is ignored so that 77 = 0, 6* = | and the system is on resonance(^ = 0). Earlier 
works on the complex sine-Gordon theory have focused only on Eqs.(5) and (6). However, 
the constraints in Eqs.(7) and (8) which are new expressions constitute an essential part 
of the SIT theory, particularly in connection with a local U(l)-gauge symmetry of SIT as 
explained later. Thus we will call Eqs.(5)-(8) as the complex sine-Gordon theory. The gauge 
symmetry structure as well as the integrability of Eqs.(5)-(8) may be best understood if we 
reformulate the complex sine-Gordon theory in terms of the action principle. The proper 
action is given by a group theoretical nonlinear sigma model action known as the deformed 
gauged WZW action defined as follows; 

S{g, A, A, (3) = Swzwig) + TT{-Adgg-^ + Ag-'dg + AgAg-^ - AA) - S^^^^^^-^^^ 

•^potential = ^ / ^""dTg-^T (9) 

where Swzwid) is the conventional SU(2)-WZW action and g is an SU(2) matrix function. 
Tr denotes the trace and T = — T = ia^ = diag(z, —i) for Pauli matrices CTj. The local gauge 
fields A = a{z, z)a3, A = a{z,z)a^ are introduced to gauge the U(l) subgroup of SU(2). 
Owing to the absence of the kinetic terms, the gauge fields A, A act as Lagrange multipliers 
which result in the constraint equations. One of the nice properties of our fomulation is that 
the equation of motion arising from the action (9) takes a zero curvature form, 

^98 = -^ + g-'dg + g-'Ag + p\T , d + A + \g''fg ]g-'5g = (10) 
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where the constant A is a spectral parameter and the square bracket denotes the commuta- 
tion. The constraint equations coming from the A, A-variations are 

5aS = ^ fTT{-dgg-' + gAg-'-A )6A = (11) 



27r 



6aS = ^ fTT{g-'dg + g-'Ag-A)6A = . 

ZTT J 



:i2) 



The action (9) is known to possess the local U(l)-vector gauge symmetry under the trans- 
form; 

g^h-^gh , A^A + h-^dh , A^A + h'^dh (13) 

where h = exp{(f){z, z)a3), as well as the global U(l)-axial vector gauge symmetry under 
g —>■ hgh for a constant h. In order to identify Eqs.(10)-(12) with the SIT equation, we fix 
the vector gauge by choosing 

A = ^T , A = (14) 

for a constant ^. Such a gauge fixing is possible due to the fiatness of A, Also, we 

parameterize the 2x2 matrix g by 

gi»7cr3gi(p(cos6'cri -sin 6*0-2 )gjJ7cr3 f ^ ^ COS If isillipe \ fl5) 

Then the parametrization in Eq.(4) arises from an identification of E,P and D with g 
through the relation 

g-'dg + ^g-'Tg-^T=(^^, , g-'fg = -z^^, _^ ) (16) 

where we have used the constraint equation (12). Also, the zero curvature equation (10) 
with the identification in Eq.(16) becomes 



tPX + i^ \ r._i_(D P 

E* -il3\~ii) ' X\P* -D 



(17) 



whose components agree precisely with the SIT equation in the sharp line limit. The con- 
straint equations (11) and (12), combined with Eq.(14), also reduce to Eqs.(7) and (8). Thus 
we have shown that the SIT equation consistently arises from the action (9) with the gauge 
fixing in Eq.(14). Moreover, the zero curvature equation (17) demonstrates the integrability 
of the SIT equation. The potential term in Eq.(9) changes into the population inversion 

r (3 r P 

"^potential = 7 ^ = 7 

which for /5 = 1 possesses degenerate vacuua at 

ip = {pn = {n -\ — )'7r, n E Z and 9 = Oq for constant . (19) 
2 
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The soliton solutions interpolating different vacuua can be obtained either by applying the 
dressing method or by using the Bacldund transformation 0. In particular, the 1-soliton 
solution is given by 

cos 09 = , ^ sech(2fe — 2hCz) 

I {a - if + 62 



Tj = {a — C,)z + {a — ^)Cz 

e = - tan"^[^^coth(262 - 2bCz)] - 2^z + 2Dz (20) 
where a, b are arbitrary constants and 

In terms of E, 

E = -2ib sech{2bz - 2bCz)e-^'^'''-^'+^''-^^^'\ (22) 

If a — ^ = 0, this solution interpolates between two different vacuua ipn and ipn+i, i-e. it 
becomes a topological 1-soliton {An = 1). On resonance where a = ^ = 0, Eq.(20) reduces 
to the 1-soliton of the sine-Gordon equation, or a 27r pulse of SIT. If a — 7^ 0, the solution 
in Eq.(20) reaches to the same vacuum asymptotically as x ±00 so that the topological 
number is zero (An = 0). Nevertheless, except for the topological number, this solution 
possesses all the properties of a soliton so that we call it a nontopological 1-soliton. It 
represents a localized pulse with a steadily varying phase. The time area of the pulse, which 
is defined by 

AS = 2 f \E\dt, (23) 



is still 271. It is important to mention that this 27r area is a mere coincidence and should not 
be confused with the 27c area of the topological one. Because of the interference between 
phases of each nontopological solitons, multi-nontopological solitons in general do not possess 
the area which is an integer multiple of 27i and the area theorem of McCall and Hahn in the 
case of inhomogenous broadening does not hold. The stability of nontopological solitons, 
unlike the topological case whose stability is due to the topological protection, arises from 
the U(l)-charge conservation law. Recall that the action (9), consequently Eqs.(5)-(8), are 
invariant under the axial vector transform g — hgh or, equivalently, 

1] rj + e for e constant. (24) 

The corresponding Noether currents and the charge are given by 

sin Lf sin Lf 



Q = / dx{J + J) (25) 



5 



where J, J satisfy the conservation law, dJ+dJ = 0. In particular, the charge of the 1-soliton 
in Eq.(20) is 

<5l-sol = -c(sign[6 ■ (a - 0] I - tan^^ — ^) • (26) 

The stability of nontopological solitons can be shown either by using conservation laws 
in terms of charge and energy as given in or by studying the behavior against small 
fluctuations 0. 

The action (9) also possesses two different types of discrete symmetries. This can be seen 
most easily in a different gauge where A = A = which is connected to the gauge in Eq.(14) 
by an appropriate vector gauge transformation. The first case is the chiral symmetry under 
the interchange, 

z ^ z and g ^ g^^{ oi i] ^ —r] , Lp ^ —Lp) (27) 

which is a characteristic of the WZW action. Unlike the case of a sigma model without the 
Wess-Zumino term, parity alone {z ^ z) is not a symmetry. In the SIT context, this is 
due to the slowly varying envelop approximation which breaks the parity invariance of the 
Maxwell-Bloch equation. The chiral symmetry generates a new solution from a known one. 
For example, the chiral transform of the 1-soliton in Eq.(20) in the resonant case = 0) is 
again a 1-soliton but with the replacement of constants a, h by 

(28) 



which changes the shape of the pulse as well as the velocity by f — > c — f . The currents and 
the charge also change into 

J ^ - J , J ^ - J , Q^-Q. (29) 

The second case is the dual symmetry of the Krammers-Wannier type|^ under 

(3 ^ -(3 , g ^ icTig . (30) 

Changing the sign of (3 makes the potential upside down so that vacuua become maxima of 
the potential and vice versa. This allows us to find a localized solution which approaches 
to a maximum of the potential asymptotically, i.e. it reaches to the completely population 
inverted state. In analogy with the nonlinear Schrodinger case, we call it a "dark soliton" . 
In practice, the dark soliton for positive P can be obtained by replacing P —(3 , z —z 
in the (bright) soliton solution of the negative /3 case. For example, the dark 1-soliton can 
be written by 



cosifie^''' = :tanh(26^ + 2bCz) 



(a - if + 62 ^(a - if + 62 

2{^a-i){z-Cz)-2iz (31) 
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Finally, we show that inhomogeneous broadening can be incorporated into our formula- 
tion naturally with minor modifications. We maintain the constraint equation (12) only and 
modify the zero curvature eqution by 



+ g-^dg + ^g-'Tg - + ~XT , B + 



(32) 



where the constant A is a spectral parameter which becomes A + ^ in the sharp line limit. 
We make the same identification as in Eq.(16) and require that g~^dg + C,g~^Tg — is 
^-independent since E is a ^-independent macroscopic quantity. This results in the SIT 
equation with inhomogenous broadening as given in Eq.(l). Once again, by using the dress- 
ing method, the 1-soliton can be obtained which is the same as in Eq.(20) but with the 
replacement 

Eq.(14) shows that each frequency ^ corresponds to a specific gauge choice in our formulation 
therefore inhomogenous broadening is equivalent to averaging over different gauge fixings. 
This implies that the inhomogenously broadened case can not be treated by a single field 
theory. Nevertheless it is remarkable that the group theoretical parametrization of E, P and 
D is still valid. Another important feature of inhomogenous broadening is that it introduces 
an anomaly term M in the U(l)-current conservation such that dJ + dJ = M and 

M = 2 cot V9[ cos{9 - 2r]) < sm{9 - 2r]) sin 2^9 > - sm{9 - 2r]) < cos{9 - 2r]) sin 2lp > 

-{cot'^ <^dr] + ^d9)d<^ ]. (34) 

This anomaly vanishes in the sharp line limit due to the constraint Eq.(8). It also vanishes 
in the case of 1-soliton and the charge remains conserved. This may be compared with the 
conserved area of topological solitons in the presence of inhomogenous broadening. It is an 
open question whether there exists a similar theorem to the area theorem concerning about 
the stability of pulses with phase variation in terms of charge and anomaly. 

We may choose different groups and coset structures for the deformed gauged WZW 
action which in our group theoretical formulation of SIT leads to other cases than the two 
level SIT. These cases, together with other aspects of SIT which were not considered in this 
Letter, will appear in a longer version of this Letter 0. 
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